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We perform a phenomenological analysis of neutrino oscillation in a four gener¬ 
ation framework introducing an additional sterile neutrino. In such a scenario, more 
than one pattern is possible that can accommodate three hieararchically different mass 
squared differences as required by the present experiments. We considered two differ¬ 
ent spectrums. Choosing the Am^s in the ranges suitable for the LSND, atmospheric 
and solar neutrino oscillation, limits on the mixing angles are derived, consistent with 
the most restrictive accelerator and reactor data as well as the atmospheric and solar 
neutrino results. The allowed mixing angles are found to be constrained very severely 
in both cases. For one mass pattern in the combined allowed zone the atmospheric 
anomaly can be explained by z/g — oscillation whereas for the other the — Vr chan¬ 
nel is preferred. The accelerator experiments CHORUS and NOMAD have different 
sensitivities in these regions and they can distinguish between the two choices. 
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1 INTRODUCTION 


The question whether neutrinos have a non-zero mass or not has remained one of the 
most tantalising issues in present day physics. In the standard model of electroweak 
theory the neutrinos are considered to be massless. But there is no compelling theoret¬ 
ical reason behind this assumption. Most extensions of the standard model allow small 
but non-zero neutrino mass. A way for probing small neutrino masses and the mixing 
between different neutrino flavours is provided by neutrino oscillations. Considering 
only two generations for simplicity, the probability that an initial z/q, of energy E gets 
converted to a after travelling a distance L in vacuum is 

= sin^26'sin^(1.27Am^L/E) (1) 


where 6 is the mixing angle in vacuum. Am^ denotes the mass difference squared in 
eV"^. LjE is in meter/MeV. The oscillatory character is embedded in the second factor 
in (|l|). The detection of this phenomenon in an experiment requires T'/L ~ Am^. 

Recently the Liquid Scintillator Neutrino Detector (LSND) collaboration has 
declared its results for a positive evidence of — z7g oscillation |^. Prior to this, 
indications of neutrino oscillations came from the well known solar neutrino problem 
and the atmospheric neutrino anomaly. 

If indeed neutrino flavour oscillation takes place, the most sensitive value of Am^ 
for detecting such a phenomenon in LSND is ~ 6eV^, which is in the right ballpark 
for the cold plus hot dark matter scenario for structure formation in the early universe 
IQ]. It remains to be seen whether the LSND results stand the test of time but already 
this has added a new impetus to the issue of neutrino mass and mixing and a number 
of investigations have been carried out recently, discussing the possible impact of this 
on particle physics, astrophysics and cosmology Q, Q]. 

The observed suppression of the solar neutrino fluxes as compared to the the¬ 
oretical predictions constitutes the long standing solar neutrino problem. A purely 
astrophysical solution to this, attributing the deficit to an inaccurate prediction of 
the fluxes by the standard solar models Q] is disfavoured by the present data 0. If 
neutrinos are massive, a plausible explanation to the solar neutrino problem is neu¬ 
trino oscillation in vacuum |Q] or the Mikheyev-Smirnov-Wolfenstein (MSW) |Q] effect 
of matter enhanced oscillations. The basic idea is flavour conversion of z/g to another 
species - active or sterile - to which the detector is not sensitive. The two generation 
oscillation explanation for the solar neutrino problem requires Am^ ~ 6 x 10~^eV^ and 
sin^ 26^ ~ 7 X 10“^ (non-adiabatic solution) and Am^ ~ 9 X 10 and sin^ 26 ~ 0.6 
(large mixing angle solution) [Q for the MSW matter oscillation to an active neutrino. 
If instead one considers oscillation to sterile neutrinos as a possible solution, the large 
angle region is excluded at 99% C.L. Q]. This region is also not consistent with the 
bound on the number of neutrino species in big-bang nucleosynthesis 0. Oscilla¬ 
tion in vacuum to an active neutrino requires Am^ ~ (0.45 — 1.2) x and 
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sm‘^2e ~ (0.6 - 1.0) O- The sterile neutrino alternative for this case is ruled out by 
the present data at 98% C.L |T2[ . 


The primary components of the cosmic-ray flux interact with the earths atmo¬ 
sphere producing pions and kaons which can decay as: 

7r+(i^=^) ^ 

^ + Uei^e) + 

These neutrinos can be detected by imaging water (Terenkov detectors, - Kamiokande 
T^, H and IMB il^ - or using iron calorimeters as is done in Frejus IM, Nusex 


TtI and Soudan2 |^. To reduce the uncertainty in the absolute flux values the usual 


practice is to present the ratio of ratios R [19 


R = 


(T )/(^^e T l^e )obsvd 


( 2 ) 


where MC denotes the Monte-Carlo simulated ratio. Kamiokande and IMB hnd R to 
be less than the expected value of unity. This deviation is known as the atmospheric 
neutrino anomaly. The preliminary results from Soudan2 agrees with this but Frejus 
and Nusex do not support this conclusion. The atmospheric anomaly, if it exists, 
can be explained by either — z/g or oscillations in a two generation picture. 

The analysis of the new multi-GeV data as well as the previous sub-GeV data of the 
Kamiokande collaboration predicts the following best-£t parameters (Am^,sin^20) = 
(1.8 xlO“^eK^, 1.0) for — oscillation and (1.6 xlO“^eK^, 1.0) for — oscillation 
IjT^. Since the required mixing angle is large, oscillation to sterile neutrinos is again 
inconsistent with the nucleosynthesis constraints 0- 


The three neutrino oscillation phenomena mentioned above - namely the so¬ 
lar neutrino problem, the atmospheric neutrino anomaly and the — Vf. oscillations 
observed by LSND group - require three hieararchically different mass ranges and a 
simultaneous explanation of all of them would need at least four generation of neutri¬ 
nos 10. LEP data reveals that there are three light active neutrino species. This is 
supported by the requirements of nucleosynthesis in the early universe. So the fourth 
neutrino has to be sterile. Introducing this new neutrino one can attempt separate 
two generation treatments for each but a more comprehensive approach would be to 
determine the parameter ranges consistent with all the experiments by a combined four 
generation analysis which can reveal the full implications of each experimental data 
on the other. In this paper we perform such an oscillation analysis in a four family 
picture. 

In a four generation scenario there are six Am^s, three of which are indepen¬ 
dent and six mixing angles neglecting the GP violating phases. We assume a minimal 
four flavour mixing scheme, in which the sterile neutrino mixes only with the electron 
neutrino, thus reducing the number of mixing angles to four. We do not make any 
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assumptions regarding these mixing angles allowing them to cover the whole range 
from 0 to 7 r/ 2 . Guided by the present data on neutrino oscillation we consider two 
different sets of hieararchical mass squared differences: 

(i) Ai3 ~ A 23 ^ A 34 = Alsnd in the LSND range, 

Ai 2 — A 24 = Aatm as preferred by the atmospheric neutrino data and 
Ai 4 = Asolar either in the MSW or in the vacuum oscillation range. 

(ii) Ai2 ~ Ai3 ~ A24 ^ A34 = Alsnd 
^23 = A atm 

Ai4 = Asolar 

Aij =1 — rrii^ \. The hieararchy in the absolute values of neutrino masses as im¬ 

plied by the above spectrums can be classified in general as 
Case (i) : (a)m^ ~ ml « m\ 


Case (ii): {a)m\ ~ ml 


ml « m\ 


or (b)m^ ~ ml » m^ ^ m^ 
or (b)mj ~ ml ~ m^ » m\ 

These are shown schematically in £g. 1. The ~ and ~ signs imply, differences by 
Asolar and A atm are neglected respectively, while the >> or << sign means dif¬ 
ference by Alsnd- Neutrino oscillation analysis remains the same for types (a) and 
(b) in both cases and one has to invoke some other experimental constraints like that 
obtained from neutrinoless double beta decay for a distinction between these [p^] . 

Models for neutrino masses and mixings assuming the existence of sterile states 
besides the three active flavours have been discussed before . After the declaration 


of the LSND results models involving extra singlet neutrinos have been constructed by 
many authors Our investigation is not motivated by any particular model. Rather, 
we do a phenomenological analysis fixing the Am^s in the ranges suitable for explaining 
the solar neutrino, atmospheric neutrino and the LSND results and determine the 
allowed values of the mixing angles, which can reconcile all the experimental data. 
The perspectives of the accelerator experiments CHORUS and NOMAD, searching for 
oscillations, are discussed in the light of our hndings. 

The plan of the paper is as follows. In section 2 we summarize the experimental 
results relevant for our purpose. In the following section the general multi-generation 
formalism for studying neutrino oscillations is developed. In Section 4 we calculate the 
survival and transition probabilities for the various experiments for the mass patterns 
(i) and (ii). The results of our analysis and some discussions are presented in section 
5. We end in section 6 with a short summary and conclusions. 


2 Experimental Results 

2.1 Laboratory experiments 

The laboratory experiments can serve as an important tool for neutrino oscillation 
search. These are either accelerator or reactor based and in general are of two types: 
(i) Disappearance experiments: in which one looks for a reduction in the initial neutrino 
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flux due to oscillation to some other flavour to which the detector is not sensitive. 

(ii) Appearance experiments: in which one searches for a new neutrino flavour, absent 
in the initial beam, which can arise from oscillation. 

Prior to the LSND results all the laboratory experiments were consistent with no 
neutrino oscillations and provided exclusion regions in the Am^ — sin^ 29 plane. In 
the region of large Am?, the sin ^(1.27Am^L/ii^) term —>■ 0.5 and = 0.5sin^26'. 
Thus the limits on sin^ 29 from the exclusion plots can be used to provide bounds 
on Puau/s- Below we summarize the laboratory experiments which give most stringent 
bounds on Pu^^p- 

The reactor experiments are suitable for a search of Pu^u^ using disappearance 
technique. The most stringent constraint is provided by Bugey which measures the 
spectrum of Vg, coming from the Pressurized Water Reactors running at the Bugey 
nuclear power plant, at 15, 40, and 95 metre using neutron detection techniques. The 
90% C.L. exclusion contour implies 1 — Pv^Ve ~ 0.05 [Q . 

Accelerators produce mainly beams - the component being small and poorly 
determined. The most restrictive accelerator based disappearance experiment mea¬ 
suring is CDHSW. The 90% C.L. bound one obtains from the exclusion plot 

presented by the CDHSW collaboration is P^ 


> 


0.95 [H 


Accelerator based appearance experiments searching for 17^ —W oscillations are 
LSND, KARMEN and E776 at BNL. The last two experiments are consistent with 
no neutrino oscillation. KARMEN has so far quoted an upper limit on the oscillation 
probability as A 3.1 x 10“^ (90% C.L.) [Q whereas from the two flavour exclusion 

areas presented by BNL one gets Pf^Fe ~ 1.5 X 10 ^ (90% C.L.) which is the more 
restrictive of the two. LSND has claimed a positive evidence for neutrino oscillations 
and reports an excess of lO.dt^gi 3.3 events over the estimated background which, 
if interpreted in terms of neutrino oscillations, corresponds to a probability Pf^Fe of 
(0.341°;20 ± 0.07)% §. 

E531 at Fermilab looks for oscillation, by Vj. appearance. Their results 

can be translated into an upper limit on the oscillation probability: Py 


~ 2 xlO-3 
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CHORUS 1^ and NOMAD at CERN search for — Vy oscillation. With 
the CERN SPS designed to deliver 2.4 xlO^® protons, CHORUS and NOMAD are 
sensitive to a minimum oscillation probability of 10“'^. 

These experiments are sensitive to different Am?, because the L/E factor is 
different for each. In Table 1 we summarize the characteristics of the experiments 
which are most constraining and of CHORUS and NOMAD whose results are awaited. 
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2.2 Atmospheric neutrinos 


Among the experiments measuring the atmospheric neutrino flux, data of most statis¬ 
tical significance have been collected by the Kamiokande and the IMB collaborations. 
For neutrinos of energy less than ~ 1 GeV, IMB finds R = 0.54 ± 0.05 ± 0.12 
agreement with the Kamiokande data R = 0.60]'^o;o5±0.05 in this energy range 


15 


m 


13, 14 


Recently the Kamiokande group has published the results of the measurement of the 
flux ratio in the multi-GeV energy range 0 - They found R = 0 . 57 lo;o 7 ± 0.07 in 
good agreement with the sub-GeV value. Another aspect of this measurement that 
can independently point towards neutrino oscillation is the dependence of R on the 
zenith-angle. The multi-GeV Kamiokande data reveals a dependence on the zenith- 
angle unlike the sub-GeV data, though the statistical significance of this result has 
been questioned 0 


For the purpose of this paper we use the sub-GeV Kamiokande 


results. 


2.3 Solar neutrinos 


At present there are four ongoing experiments that are measuring the flux of solar 
neutrinos; the Homestake ^'^Cl neutrino capture experiment, the Kamiokande u-e scat¬ 
tering experiment and the neutrino capture experiments of the SAGE as well as 
the GALLEX collaborations. The measured rate for for the chlorine experiment 
is 2.55 ± 0.17 (stat) ± 0.18 (syst) SNU as compared to the standard model |0| pre¬ 
diction : 9 . 3 li ;4 SNU including metal and helium diffusion and 7.0l5;o SNU without 
any diffusion. The observed rates in the Ga experiments GALLEX 0] and SAGE 
are respectively 79 ± 10 (stat) ± 6 (syst) SNU and 69 ± 11 (stat) I 7 (syst) SNU 
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while the theoretical prediction is 137 I 7 SNU including diffusion and 126^6 SNU for 
no diffusion. The measured flux in Kamiokande |0] is [3.0 ± 0.41 (stat) ± 0.35 (syst)] 
xl0“® This is 0.45 of the theoretical predictions including diffusion and 0.61 

of the standard model rate without diffusion 0]. Thus all the experiments indicate 
deviations from the standard model predictions. The degree of depletion differ from 
experiment to experiment. Since each type of experiment is sensitive to different parts 
of the solar neutrino energy spectrum, it is plausible that the suppression mechanism 
is energy dependent. It has remained an unsettled issue for a long time and new 
experiments [0 are being set up to examine it from various angles. 


3 Neutrino Oscillation 

3.1 Vacuum Oscillation 

If the neutrinos are massive then quark mixing suggests a leptonic mixing matrix 
analogous to the GKM matrix. Then, the flavour eigenstates, i.e. the states produced 
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in weak interaction decays are not the same as the mass eigenstates (which are the 
neutrino states that propagate) but linear combinations of them. 

z/„ = (3) 


U is the unitary mixing matrix. This leads to the possibility of neutrino oscillations. 
For N neutrino generations the expression ( 0 ) generalises to 


Pu. 


tUp — ■ 


tiL 

4 UaiUpiUajUpj sin (— ) 


(4) 


i, j varies from 1 to N for N generations. = 2A7m{Ey /MeV){eV‘^/ /S.ij). The actual 

form of the various survival and transition probabilities will depend on the spectrum 
of Am^ chosen and the explicit form of U. 

If Am^ is such that a particular A >> L, then the corresponding oscillatory 
term sin^ ttL/A ^ 0, whereas A << L would imply a large number of oscillations and 
consequently the sin^7rL/A term averages out to 1/2. 

Since CP violating phases are discarded, U is real and in the minimal mixing 
scheme, is a function of four angles. One can express it in general as the product 
of four 4x4 rotation matrices, Rij. The order in which the multiplication is to be 
performed is arbitrary leading to several different forms of U. In the quark sector 
considerable effort has been devoted to specify what would be a ‘good choice’ for N=3 

All choices being mathematically equivalent, a good choice 


as well as N > 3 38 


is mostly a matter of convenience, depending on how directly one can connect the 
experimentally determined quantitities with the mixing angles. In the neutrino sector 
the experiments measure the energy and detector cross-section averaged survival or 
transition probabilities. Consequently a judicious choice of U in this case would be 
that in which one can minimise the number of mixing angles appearing in the various 
survival and transition probabilities. With this as the guideline for determining the 
mixing matrix U, we observe that the order in which the product of the rotation 
matrices are to be taken depends on the Aijs chosen. As a general principle we £x 
the following rule: the rotations are to be performed in the order of increasing mass 
hieararchies. This is true even for a three generation analysis. As a consequence of 
CPT and CP invariance there is no distinction between Pu^up, PvaVp , Pvpvc > Pupv^ 
our analysis. 


3.2 Matter Oscillation 

The discussion on matter oscillations is geared towards the solar neutrino problem 
since we neglect the matter effects for the atmospheric neutrinos. For the mass scales 
involved in our analysis this is a good approximation [^ . 

We have kept A14 in the range suitable for solving the solar neutrino problem for 
both the mass patterns, implying z/g — Vg oscillation as the dominant mode for depletion 
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of solar neutrinos. In a combined analysis, mixing of v^. with and is also expected 
to affect the probabilities. 

Solving the neutrino propagation equations in matter for more than two genera¬ 
tions and arbitrary values of neutrino masses is in general a non-trivial exercise. MSW 
analysis for three neutrino generations and the conditions under which it simplihes 
have been studied by many authors |^. A particularly simplifying assumption is one 
in which the problem reduces to an effective two generation case ^ . Oscillation to 

sterile neutrinos as a solution to the solar neutrino problem has been pursued by many 
authors in a two generation framework - for both vacuum , and matter oscillations 


4^ . In the presence of more than two generations the propagation equations in matter 
differ from the two generation Vf. — case, due to the charged current interactions of 
only Vf. with matter as well as the different neutral current interactions of the three 
active flavours and Vg. In our analysis we follow the notations of 0 but extend it 
for four generations including the changes in the effective potentials arising from the 
presence of sterile neutrinos. 

The neutrino propagation equation in matter for four neutrino generations is. 


dx 


( ^e\ 
Vi, 
Vr 

\VxJ 


—Ml 
2E ^ 


( Ve\ 
Vf, 
Vr 

\Vx J 


(5) 


where, the mass matrix in flavour basis is 


Mn^ = 


Mf 

= UMo^Ui 

+ Ma 

/o 

0 


0 

0 

\ 

0 

A21 


0 

0 


0 

0 


A31 

0 


VO 

0 


0 

Ai4 / 


/Ai 

0 

0 




0 

0 

0 

0 



0 

0 

0 

0 



1 0 

0 

0 

A2 ) 





Aol 


( 6 ) 


(7) 


( 8 ) 


Ai = 2\/2GFnf.E and arises from the charged current interaction of Vf, with electrons 
while A 2 = V^GpnnE is due to the neutral current interaction of z^e, as well as Vr 
with the neutrons. Gp is the Fermi coupling constant. Ug and are the ambient elec¬ 
tron and neutron density respectively. E is the neutrino energy. For A12, A13 >> Ai, 
can be separated into a 2x2 neutrino submatrix exactly of the two generation form 
and two decoupled neutrinos. Under these conditions MSW resonance occurs between 
the hrst and the fourth generation while the second and the third generation remain 
unaffected by matter. This conclusion has been checked in [M for one mass squared 
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difference in the atmospheric range. Since Alsnd >> Aatm, the above arguement is 
expected to be valid for this mass scale also. We evaluate the survival probabilities in 
this approximation in the subsequent section for both the mass patterns. 


4 Survival and Transition Probabilities 


The oscillations are characterised by three oscillation wavelengths Xlsnd, Xatm, Xsolar 
corresponding to the three mass scales, Alsnd, ^atm and Asolar in the problem. 
We note the following general points: 

(i) Accelerator and Reactor Neutrinos 

For the accelerator and reactor neutrinos, the energy and length scales are such that 
XsoLAR and Xatm » L (see Table 1) and the oscillations driven by these mass scales 
are absent. Thus the one mass scale dominance often used in the context of accelerator 


and reactor neutrino oscillation ^ is a valid approximation, the oscillations being 
driven by Xlsnd- We further note that for Bugey Xlsnd « L so that sin^{ tiL/Xlsnd) 
averages to 1/2. 

(ii) Atmospheric neutrinos 

For the atmospheric neutrinos in the energy range ~ (0.1 - 1) GeV travelling through 
a distance ranging from ~ (10 - 10^) km, Xlsnd « L and siv?{TiL/X lsnd) can be 
replaced by the average value 1/2. On the other hand Xsolar » L and the oscilla¬ 
tions driven by the mass scale relevant for the solar neutrino problem vanish in this 


case. 


(iii) Solar neutrinos 

For oscillation of solar neutrinos in vacuum Xlsnd and Xatm « L and the terms 
involving these average out to 1/2. 


4.1 Mass Spectrum (i) 

Following the prescription described in section 3.1, 


U 


R 23 R 13 R 12 R 14 

/ C12C13C14 

—C23S12C14 — S13C12S23C14 
'Sl2'S23Cl4 — 'S13C12C23C14 

V -Sl4 


SI2C13 Sl3 

C12C23 — 'Sl3'Sl2'S23 'S23C13 

— Cl2'S23 — 'S13'S12C23 C13C23 

0 0 


C13C12SI4 \ 

—C23S12S14 — 513012523514 
512523514 “ 513C12C23S14 

Cl4 / 

(9) 


where = cos 6ij and = sin6*jj, here and everywhere else in the paper. Let us now 
see what the probabilities for the various experiments are in this case. 

(i) Accelerator and Reactor Experiments 


ft.F. = 1 - 2sl,cl, (Bugey) 


(10) 









= 1 - 45^34(1 - cl,sl,)sm\nL/\LSND) {CDHSW) (11) 

Pv^v,= ^c\^s\^s1^siv?{'kL/\lsnd) {LSND,E776) (12) 

= 4444sin2(7rL/AL5^fl) {E531, CHORUS, NOMAD) (13) 

(ii) Atmospheric neutrinos 

The general expression of (@) for N flavours in terms of the neutrino transition and 
survival probabilities is 


R 


Ill'll T "^McPuiil/e 

F + ^ F 


(14) 


where tmc = ipe + as obtained from a Monte-Carlo simulation. In this 

case the relevant probabilities appearing in ( 0 ) are 

Pu,u, = 1 - 24si3 - AcI^sI^cI^AyiS^tiL/Xatm) (15) 

= 244'S23+4c?3Cl2Sl2(c23Ci2-Si3Si2S23)(c23Sl2 + Sl3Ci2S23)sin^(7rL/AArM) (16) 

Pyii^ii — 1 ~ 2 C 33 S 23-|-244 ~ 4sin^(7rL/AArM)(c23Cl2 — Sl3Sl2-S23)^(c23'Sl2 + Sl3Cl2'S23)^ 

(17) 


(iii) Solar neutrinos 


(a) Vacuum Oscillations 

The electron neutrino survival probability in this case is, 

-Pi/ei/e = O3C12-P2VAC + S13 -f S12C13 ( 18 ) 

where P 2 VAC is of the form of the two generation vacuum oscillation probability: 

P 2 VAC = 1 - sin‘^29u sin'^ (tiL/Xsolar) (19) 


(b) MSW oscillations 

If the mass hieararchies, mixing angles and the density distributions are such that 
one has resonance between the hrst and the fourth mass eigenstates, while the second 
and the third mass eigenstates remain independent of matter density then the mixing 
matrix in matter is Um = R 22 ,RizRi 2 Riam■ The mixing between the hrst and the 
fourth generation gets modihed by the matter effects as. 


t(in26i4rn — 


Ai 4 sin 2014 

Ai4 cos 2014 — 2y/2GpneffE 


( 20 ) 


where we dehne Ug// as 

Ueff = Ci 3 Ci 2 Ue — — 77 ,,^ 


( 21 ) 
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1 this reduces to the two generation limit for oscillation between 
From (1^), the resonance condition between the hrst and the fourth 


In the limit C 13 , C 12 
Vf. and Vg 

generation in the presence of the other two generations becomes 

2y/2GpTi^fj^E = /\i4^cos29i4^ 


( 22 ) 


The difference of this with the three generation resonance is to be noted. While the 
three generation resonance condition gets modihed by one additional mixing angle 
here the mixing angles with the second as well as the third generation appear in the 
resonance condition through rie//. The calculation of the survival probability is 
then a straight-forward generalisation of the standard two or three generation MSW 
scenario and one gets, 

T’l/gi/g = Ci2Ci^PmSW + ■S 13 -f- <S]^2 Ci3 (23) 

Pmsw ~ 0-5 + (0.5 — 6[E — Ej^')X')cos29i4^cos29i4^]\/[ (24) 

Ea = Ai4Cos29i4,/2^/2GFneff gives the minimum u energy that can encounter a reso¬ 
nance inside the sun. For this case the jump probability X, between the hrst and the 
fourth mass eigenstates, as obtained in the Landau-Zener approximation is 

1 (272 f f 

X = exp (— 7 rAi 4 sin^ 26 *i4/4F^cos26'i4 -(— p^)res) (25) 

rieff ax 


We note that unlike the case discussed in [^, the mixing angles C 42 , c ^3 appear in the 
expression of the jump probability, via rie// dehned in (pi]). 


4.2 Mass Spectrum (ii) 


This case has been considered also in following slightly different notations. We 


discuss this to make the paper self contained and to facilitate the comparison between 
the results obtained for the two mass patterns. Following the rule we specihed U in 
this case is 


U — R^^R^oRo'tR 


13-rri2-n.23-ri-i4 

/ C12C13C14 S12C13C23 — S13S23 C13S12S23 + S13C23 C13C12S14 ^ 

—S12C14 C12C23 C12S23 —S12S14 

S13C12-S14 —S13S12C23 — C13S23 —S12S13S23 + C13C23 —S13C12S14 

—S 44 0 0 C 44 


(26) 


This reduces to the three generation case discussed in in the limit of 6*14 going to 
zero. We next calculate the probabilities explicitly for the various experiments. 

(i) Accelerator and Reactor Experiments 


P- - — 1 — 2r^ P 


+ 2 C 13 C 42 


{Bugey) 


(27) 
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Pv 


= 1 - sin^ 2ei2sm^ {'kL/Xlsnd) (CDHSW) 


Pu^u. = ^c{^s{2ci^sm.\'KL/\LSND) iLSND,E776) 
p--^ = Ac\^sl^sl^sm^{TiL/\LSND) (^531, CHORUS, NOMAD) 


(28) 

(29) 

(30) 


(ii) Atmospheric neutrinos 


For the chosen mass pattern and mixing the probabilities appearing in (|T^ are 

^2 , 0 . 4^4 4(ci3Si 2C23 - Si3S23)^(ci3Si2S23 + Si3C23)^sin^(7rL/AATM) 


— r ZC]_3C]_2 


n3'-i2 

P — 20 0 0 

j - ZiCiqC-ioO 


— ^'-13'-12'^12~4C]^2C23'S23(ci3’Si 2C23—Si3S23)(ci3Si2S23 + Si3C23)sin (ttL/Xatm) (32) 


Pu..u.. = 1 


2ci2'Si 2 4c]^2C23'S23sin [ttL/ Xatm) 


(31) 

(32) 

(33) 


(iii) Solar neutrinos 


(a) Vacuum Oscillations 

The electron neutrino survival probability in this case is, 

PueUe = Ci^Ci2P2VAC + (S 12 C 13 C 23 — 513823 )^ + (C 13 S 12 S 23 + 513023 )^ (34) 


where P 2 VAC is given by eqn. (pTf ). 

(b) MSW oscillations 

In this case the mixing matrix in matter is, Um = -Ri 3 -Ri 2 -R 23 -Ri 4 M- For this also the 
resonance can be assumed to happen between the first and the fourth mass eigenstates 
whence the mixing angle 0 i 4 m and the resonance condition continue to be given by 
eqns. ([^) and (^^. The probability however is different and can be expressed as. 


Pue^e — Ci2Ci3-Fm5W + (S12C13C23 — S13S23) + (C13S12S23 + S13C23) 


(35) 


where Pmsw is defined in eqn. 


5 Results and Discussions 


For both the mass spectrums discussed in sections 4.1 and 4.2, the survival or tran¬ 
sition probabilities for the accelerator and reactor experiments are functions of any 
two of the mixing angles - 6 ^ 12 , 6 ^ 13 , 6*23 - and one mass squared difference, Alsnd- In a 
realistic analysis one has to average the probabilities over the LjE distributions of var¬ 
ious experiments and fold it with the detector cross-sections. We adopt the approach 


followed in since one mass scale dominance is a good approximation in our case 
also. As noted in [45| in this limit one can make a one to one correspondence between 
sin^ 26 as obtained from a two flavour analysis and the angular factor appearing in a 
three generation calculation. If we now £x AmS around Alsnd ~ 6 eV^ then from the 
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bound on sin^ 29 at this Am^ from two flavour exclusion contours one can constrain the 
three generation mixings. While the other experiments present the exclusion contours 
at 90% C.L., LSND gives their plots at 95% C.L. |^. Thus for this we use the quoted 
value for the probability at 90% C.L. and use the limit sin^( 7 rL/AL 5 Arz)) ^ 1. 

The probabilities in the atmospheric neutrino case are functions of the three mix¬ 
ing angles 6 ^ 12 , 6 * 13 , 623 and one mass squared difference. In our analysis we approximate 
the {tiL/Xatm) factor by its averaged value 0.5 as is often done in the context of 
the sub-GeV data Q. This can be improved by an averaging over the incident 

neutrino energy spectrum, the zenith-angle of the beam as well as the final lepton en¬ 
ergy tmc is taken to be 0.45 from a detailed Monte-Carlo simulation including 


the effects of muon polarisation [^. Using the sub-GeV Kamiokande results: 


0.48 <R< 0.73 (90%C.L.) 


(36) 


We determine how much of the area, admissible from neutrino oscillation searches at 
laboratory, permits R to lie in the above range for hxed values of the third mixing 
angle. The allowed ranges of this mixing angle consistent with the atmospheric data 
are also obtained. 

Finally we check the compatibility of the mixing angles found from this combined 
analysis with the solar neutrino results. The sensitivity of CHORUS and NOMAD in 
these areas are explored. Below we discuss the results for each mass spectrum. 


Mass Spectrum (i) 

The various constraints in this case from the accelerator and reactor results are: 


2 s? 3 Ci 3 < 0.05 (Bugey) 

(37) 

444(1 - 44 ) < 0.08 (CDHSW) 

(38) 

45^34^23 < 0.003 (U776) 

(39) 

0.0002 < 4s?34s^3 < 0.0069 (LSND) 

(40) 

4.S23G3'^i3 ^ 0.02 (77531) 

(41) 

From the Bugey constraint (|37|) 


s ?3 ~ 0.026 or s ?3 ~ 0.974. 



Conhnining within these limits we scan the whole range 0 < A 1 to determine 
the admitted area in the 3^3 — plane consistent with the accelerator and reactor 
experiments listed above. A large portion of the parameter space can be ruled out 
combining all the restrictions and three allowed sectors are obtained - presented in 
hgs. 2 a, 2 b and 2 c respectively. 
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Fig. 2a shows the zone where 0.974 ~ sfg < 1.0. As is clear from the hgure, the 
most stringent limit on S 23 , in this range, comes from E776 and LSND constraints. The 
permitted area of £g. 2a corresponds to the case —> 1. In this regime the relevant 
probabilities for atmospheric neutrinos given by eqns. (0 to ( 0 ) can be expressed 
as ~ 1, ~ 0, Py^ ^ ^ ^iu^ {29 12 T 2023). This is the oscillation 

limit, where consistent solutions to the atmospheric puzzle may be obtained. However 
for the solar neutrino survival probability, eqn. (|^) , the limit 1 would imply 

that the coefficient of the vital term responsible for the MSW effect, Pmsw, becomes 
very small. Consequently Py^y^ —>■ 1 due to the factor sfg in eqn. (^) contrary to 
the results from the solar neutrino experiments discussed in section 2.3. Similarly for 
the vaccuum oscillation probabilitiy (eqn. (|T^)) also one would require that the factor 
c^gC^g multiplying the energy and Am^ dependent term, P 2 VAC should not be too small 
which would again prefer low sfg. Thus this area is disfavoured by the solar neutrino 
data. 


For low s^g 




0.026) the allowed values of S 23 is severely constrained by E531 


result and only very high (> 0.99) or very low (< 0.01) S 23 are admissible. 


Fig. 2b shows the low s^g - high S 23 zone. In this regime E776 data puts a 
stronger constraint on s^g than Bugey and sets the limit s^g < 0.001. In this part of 
the parameter space a consistent solution to the atmospheric neutrino anomaly cannot 


be found. This region corresponds to sjg —>• 0, s 


23 ^ J-, whence from eqns. (|T^ to 

oscillations 


0) Pu,u, ^ 1 - 2 cf 2 S? 2 , - 0 and Py^y^ ~ 1. Thus in this zone Ue - v-. 

take place driving the ratio of ratios in a direction opposite to that required. 

Fig. 2c is the low s^g - low S 23 region. In this portion of the parameter space most 
severe restrictions are from E531, LSND and Bugey giving the following constraints: 
0.01 < sfg < 0.026, 0.002 < S 23 < 0.005. The angles in this sector can be approximated 
as Si 3 —0, S 23 —> 0. In this limit the relevant probabilities for the atmospheric 
neutrinos (eqns. ([T5|) to (|^) are Py^y^ ~ 1 - 2 c^ 2 'Si 2 > Pv^v^, - 2 cf 2 sf 2 and Py^y^ ~ 
1 — 2 c^ 2 ^i 2 - Thus this belongs approximately to the oscillation region where the 

atmospheric anomaly can be explained consistently. Varying sfg and S 23 in the range 
shown in Eg. 2c the limits on 5^2 compatible with the atmospheric neutrino constraint 

(|3^) are presented in £g. 3. Two allowed bands are obtained for 5 ^ 2 : 0.10 ~ 5^2 ~ 0.37 
and 0.61 ~ s ^2 ~ 0.88. Taking six different s^ 2 ®) from the above range corresponding 
to the edges of the two permitted zones and the parts inside them, the allowed area in 
the sfg — S 23 plane from accelerator, reactor and atmospheric neutrino data is shown 
shaded in £g. 4. Compatibility with the solar neutrino data disfavours the high values 
of 8^2 and only part of the parameter spaces shown in hgs. 2c, 3 and 4 satisfy all the 


experimental constraints. 

Fixing ~ 6 eV^, from eqn. ([I^) Py^y^ = O.lOs^gC^gC^g for CHORUS and 

NOMAD. For typical values of mixing angles from this combined allowed zone, namely. 


’13 


= 0.003 and s^g = 0.02, P, 




3.12 X 10 which is greater than the minimum 
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sensitivity, 10 of these experiments and can be marginally within their reach. 


Mass Spectrum (ii) 


In this case the constraints on mixing angles at Alsnd ~ 6 eV"^ 
laboratory experiments are: 

from the various 

2 c? 3 Ci 2 - 2 C 13 C 12 < 0.05 {Bugey) 

(42) 

4c^2S?2 < 0-08 {CDHSW) 

(43) 

4012^12^13 < 0-003 (P776) 

(44) 

0.0002 < 4c^i2Si2Ci3 < 0.0069 {LSND) 

(45) 

4oi2^12^13 < 0.02 (P531, from —s- v^) 

(46) 


For this case also three allowed areas in the relevant 3^2 — plane are obtained, from 
the laboratory constraints. The experiments which are most restrictive in this case 
can be different in general from the mass spectrum (i). These regions are displayed in 
hgs. 5a, b and c respectively. The CDHSW constraint (^31) gives, 

~ 0.02 or ~ 0.98. 


Fig. 5a shows the area for which 5^2 ~ 0.98. Allowed ranges of in this region 
are determined by E776 and LSND. This corresponds to 5^2 ^ ^^is limit for the 
atmospheric neutrinos — 0 and — 1, from the expressions (|32[) and (^). 
Thus this is the — Vt oscillation regime and is not consistent with the atmospheric 
anomaly. 


For s 
> 


< 

rsj 


0.02, the Bugey constraint 
to be ~ 0.97 or ~ 0.03. 

Fig. 5b contains the area where SI 2 ~ 


restricts the permissible values of Sgg 


< 


0.02 and s 


*12 


13 


> 


0.97. In this region s 


13 


0. Then, the probabilities given by eqns. 


neutrinos assume the following forms: Pi 




— ^'- 23*235 


to (1^ for the atmospheric 

— 2C2gS2g and Py„v„ ~ 


1 — 2c|gS2g. Thus, in this limit Pi/ — 0 and the atmospheric puzzle can be explained 


by Vf. — Vfj, transitions. However, since sfg ^ 1, this region is incompatible with the 
solar neutrino flux measurements as is evident from (|3^ and (^). 

Fig. 5c shows the area 0 < Sg 2 < 8 x 10“^ and 0 < < 0.03 


’13 


.2 

’13 


0 , 


’12 


Here, E776 data 
is determined by 
0. In this range eqns. 


puts a tighter bound on 5^2 than CDHSW. The constraint on Sgg 
Bugey. This region corresponds to the limit s 
to ( p3D for atmospheric neutrinos imply 
So this corresponds to —z/,- oscillations in the atmosphere. Substituting these in (|36D 
one gets the limits on s|g as 0.162 < s^g < 0.838. A similar situation was discussed in 


1 P ~ (1 and P 


^2 .2 
- 23 '^ 23 - 


47| , which considered the three generation limit of the mixing matrix 
for all §23 lyiiig within this limit the whole of the parameter space shown in Eg. 


As in 1^ 


5c 
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is consistent with the condition (^). Thus in hg. 6 we present the allowed region in 
the 5^2 ^ plane consistent with accelerator, reactor and atmospheric neutrino data 


’13 


choosing one representetive value of S 23 (= 0.2) from the above range. The difference 


with 


is, here, E776 constraints further narrows down the allowed range of s^ 2 - 
Since in this zone both 5^2 close to zero this is consistent with the solar 

neutrino results. 

For CHORUS and NOMAD from eqn. (l30|) = 0.16s^2Ci2'Si3- Taking two 

typical values of 3^2 and s ^3 , namely, s ^2 = 10”"^ ^'^cl 5^3 = 0 . 02 , from the combined 
allowed zone, Pu^u-r ~ 10 “® which is below the minimum sensitivity attainable in these 
experiments. 

We observe that in both cases, the solar neutrino survival probailities eqns.(|T^ 
and for the vacuum oscillation case and eqns. (^) and (^5]) for matter oscillations 
depend on A 14 , ^14 as well as on combinations of ^ 12 , 6 * 13 , 623 which depend on the mass 
pattern. Fixing the values of these mixing angles in the region determined by the 
atmospheric and laboratory results one can End the allowed area in the A 14 — sin^ 26*14 
plane using the solar neutrino data. The solar neutrino probabilities in this case differ 
from the two generation z/g — I's oscillations case, due to the presence of the mixing 
angles with the other generations. Also, in the two generation case involving just z/g 
and Ug Kamiokande would be sensitive to z/gS only but here one has the additional 
possibility of a simultaneous transition to z/^s as well as which can interact in the 
Kamiokande detector by virtue of their neutral current interactions. It has been shown 
in that the appearance of one mixing angle in the expression of probability 

leads to a larger area in the MSW parameter space. Similar conclusions might be 


obtained here also. The following points of differences are to be noted: 

113 or m three generations are involved and the resonance condition as well 


In 


as the survival probabilities are affected by the presence of only one additional mixing 
angle, whereas here one has two or three such mixings. 

(ii) In the above references oscillation between active species were considered and hence 
the term involving the neutron density was absent and the jump-probability retained 
its two generation form. Here due to the asymmetric interaction between the active 
and the sterile species the jump-probability between the first and the fourth state is 
also affected by mixing with the other generations. 

A definitive prediction regarding how these would change the two flavour allowed 
zones need a detailed numerical analysis of the solar neutrino data including the neu¬ 
trino fluxes, density profile, interaction cross-sections and a thorough treatment of the 
theory errors and their correlations. This is not performed here. 

Similarly in the vacuum oscillation case also one can probe whether the presence 
of the other mixing angles will alter the two flavour parameter space. 

We note that for the mass spectrum (ii) in the combined allowed zone depicted 
in fig. 6, both 5^2 and S 43 stay close to 0. Thus in this case, from eqns. (|3^ and 
1 ) for the vacuum and matter oscillation case respectively, the presence of the other 
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mixing angles is not expected to change the two flavour parameter space signiflcantly. 


6 Summary and Conclusions 

We have performed a combined analysis of the accelerator, reactor, atmospheric and 
solar neutrino data in a four generation framework introducing a sterile neutrino, Ug. 
In such a scenario there are in general six mass squared differences, three of which 
are independent and six mixing angles, neglecting CP violation in the lepton sector. 
We assume that Vs mixes only with z/g, thus reducing the number of mixing angles 
to four - 6 ^ 12 , 6 * 13 , 6^23 and 9u. Fixing the three independent Am^s around the ranges 
from two generation analyses of the LSND, atmospheric and solar neutrino data, we 
determine the mixing angles consistent with all the experimental constraints. We 
consider a picture where A 14 is fixed in the solar neutrino range (either MSW or vacuum 
oscillation). Then one can think of two different mass patterns for the remaining five 
Am^s - the mass spectrum (i) in which two Am^s are in the atmospheric range and 
the other three in the LSND range and the mass spectrum (ii) where one Am^ is in 
the atmospheric range and the remaining four in the LSND range. For both cases 
one can parametrise the mixing matrix in such a way that the probabilities for the 
accelerator and reactor experiments are functions of only two mixing angles and one 
independent Am^ viz Alsnd - Fixing Alsnd at 6eD^ we map out the allowed zone in 
the sin^ 6*13 ~ sin^ 6*23 (sin^ 6*12 — sin^ 6 *i 3 ) plane for the mass spectrum (i) ((ii)). Using 
the atmospheric neutrino constraint the above area can be further restricted and the 
permissible ranges for the remaining mixing angle can be determined. Next we examine 
whether the combined allowed area thus obtained is compatible with the solar neutrino 
results. 

In general for both mass patterns the following picture emerges: the accelerator 
and reactor experiments give three allowed sectors of relevant mixing angles. In two 
of these zones a simultaneous solution to the atmospheric anomaly is possible. One 
among these is disfavoured by the solar neutrino data - leaving us with a narrow range 
for permitted ^12, 6*13 and 6*23. For the mass spectrum (i) the admitted zone from all the 
input information is the one where the atmospheric puzzle can be explained by z/g — z/^ 
oscillation, while for the mass spectrum (ii) in the combined allowed zone it is due to 
— W oscillations. For both mass patterns Pui,u-r foi' CHORUS and NOMAD is much 
below the minimum reach of these experiments in the region where the solution to 
the atmospheric problem is via z/^ — z/,- oscillations. For mass spectrum (ii) this being 
the combined allowed zone, cannot be explored by CHORUS and NOMAD whereas 
for (i) in the favoured zone for CHORUS and NOMAD can be greater than the 
minimum sensitivity of 10“^ and could be probed by these experiments. 

In conclusion we would like to mention that though a three flavour mixing scheme 
cannot accommodate the three hieararchically different mass ranges required for LSND, 
atmospheric and solar neutrino oscillations, in a four generation framework with an 


16 



additional sterile neutrino there are more than one possible mass spectrums that can 
account for all the data simultaneously. We discussed two such mass spectrums. In 
both cases our analysis assumes that the solar neutrino oscillation is driven mainly by 
Ve — ^s transitions. If this scenario is conhrmed by the future solar neutrino experiments 
and the other experimental inputs do not change signihcantly as more data accumulates 
then it is necessary to go to a four generation picture. We have shown that the 
implications of CHORUS and NOMAD are different in the two cases and thus they 
can distinguish between these. 

In this article the allowed areas are obtained by hxing Alsnd ~ 6 eU^. LSND is 
sensitive to the range 1-10 eU^ and it remains to be seen what best-£t value, consistent 
with KARMEN and BNL-E776, emerges when more data is accrued. We beleive that 
the general conclusions obtained in this analysis will remain the same for any other 
value of Am^ in the above range, though the precise values of the allowed mixing angles 
may be different. 
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cussions, a careful scrutiny of the manuscript and encouragement at every stage of this 
work. She also wishes to thank Dr. Kamales Kar for discussions, help and encourage¬ 
ment. Financial support from the Council of Scientific and Industrial Research, India 
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Table 1: The characteristics of the most restrictive accelerator and reactor experi¬ 
ments. Xlsnd and Xatm are calculated for A^snd ~ 6 eV"^ and Aatm ~ 10“^ 
respectively. 


Experiment 

E 

L 

Xlsnd 

Xatm 

Bugey 

~ 5 MeV 

~ 40 m 

~ 2.08 m 

1250 m 

CDHSW 

2<E <20 GeV 

~ 1 km 

(0.83 - 8.33) km 

(500 - 5000) km 

E776 

1 - 10 GeV 

~ 1 km 

(0.416 — 4.16) km 

(250 - 2500) km 

E531 

~ 50 GeV 

0.949 km 

~ 22 km 

~ 12500 km 

LSND 

(36 - 60) MeV 

30 m 

(15-25) m 

(9 - lh)km 

CHORUS/NOMAD 

30 GeV 

0.8 km 

12.5 km 

7500 km 
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FIGURE CAPTIONS 


Figure 1: The level diagrams showing the possible mass hieararchies (not to scale). 
Figure 2: The allowed region in the sfg — S 23 plane from accelerator and reactor data 
for the mass pattern (i) 

(a) The region between the solid lines is allowed by LSND; the area to the right of the 
big-dashed line is permitted by E776 and that below the small-dashed line is allowed 
from CDHSW. 

(b) The admiited area from LSND is between the solid lines while that from CDHSW 
and E-531 is above the small-dashed and medium-dashed lines respectively. The area 
to the left of the big-dashed line is allowed from E776. 

(c) The area between the solid lines is allowed by LSND while that below the big- 
dashed, small-dashed and medium-dashed lines are allowed from E776, CDHSW and 
E531 respectively. 

In each of these the combined allowed area is marked ‘allowed’. 

Figure. 3: The allowed region in the <s ^2 plane consistent with accelerator, reactor 
and atmospheric neutrino constraints. S 33 is varied in the range determined from £g. 
2 (c). 

Figure 4: The allowed area of fig. 2c that is consistent with the atmospheric neutrino 
constraint is shown shaded for six different values of S 32 from the admiited range in 
fig- 3. 

Figure 5: The allowed region in the s ^2 ~ plane from accelerator and reactor data 
for the mass pattern (ii) 

(a) The area between the curved solid lines is allowed by LSND, while that to the right 
of the vertical solid line is consistent with the CDHSW constraint; the zones above the 
small-dashed, medium-dashed and big-dashed lines are allowed from Bugey, E531 and 
E776 respectively. 

(b) The permiisible area from LSND is between the curved solid lines; the areas above 
the small-dashed, medium-dashed and big-dashed lines are allowed from Bugey, E531 
and E776 respectively; the region to the left of the vertical solid line is admitted from 
CDHSW. 

In each of these figures the area marked as ‘allowed ’ is consistent with all the con¬ 
straints. 

Figure 6 : The allowed region of fig. 5c that is consistent with the atmospheric neutrino 
constraint is shown shaded. 


19 



References 


[1] C. Athanassopoulos et al, |nucl-ex/9504002| (1995); W.C. Louis, Nucl. Phys. (Proc. 
Suppl.) B38, 229 (1995). 


[2] J.R. Primack et al, Phys. Rev. Lett. 74, 2160 (1995). 


[3] G. Raffelt and J. Silk, |hep-th/9502306 . 

[4] See for example G.M. Fuller, J.R. Primack and Y.Z. Qian, astro-ph/950208i ; D.O. 
Galdwell and R.N. Moliapatra, Preprint No. UGSB-HEP-95-1, hep-ph/9503316| ; 
S.M. Bilenky et. a/.. Preprint No. DFTT 25/95, JHU-TIPAG 95013, hep-ph 
95044051 (1995). 


[5] J. N. Bahcall and M. H. Pinsonneault, Rev. Mod. Phys. 64, 885 (1992); S. Turck- 
Ghieze and 1. Lopes, Astrophys. J. 408, 347 (1993); See also J. N. Bahcall and 
M. H. Pinsonneault, Preprint No. lASSNS-AST 95/24, to appear in Rev. Mod. 
Phys. October 1995. 


[6] N. Hata and P. Langacker, Phys. Rev. D47, 2220 (1993); J.N. Bahcall, Phys. Lett. 
B338, 276 (1994); W. Kwong and S.P. Rosen, Phys. Rev. Lett. 73, 369 (1994); 
S. Parke, Phys. Rev. Lett. 74, 839 (1994). 

[7] For a review on vacuum oscillation see S.M. Bilenky and B.M. Pontecorvo, Sov. 
Phys. Usp. 20, 776 (1977). 

[8] L. Wolfenstein Phys. Rev. D34, 969 (1986); S. P. Mikheyev and A. Yu. Smirnov, 
Sov. J. Nucl Phys. 42(6) 913 (1985); Nuovo Gimento 9c 17 (1986). 

[9] N. Hata and P. Langacker, Phys. Rev. D50, 632 (1994) 

[10] X. Shi, D. Schramm and B. Fields, Phys. Rev. D48, 2563 (1993). 

[11] N. Hata, Univ. of Pennsylvania Preprint No. UPR-0605T, 1994. 

[12] P.I. Krastev and S.T. Petcov, Phys. Rev. Lett. 72, 1960 (1994). 

[13] K. S. Hirata et al, Phys. Lett. B280, 146 (1992). 

[14] Y. Fukuda et al.., Phys. Lett. B335, 237 (1994). 

[15] D. Gasper et a/., Phys. Rev. Lett. 66, 2561 (1991); R. Becker-Szendy et a/., Phys. 
Rev. D46, 3720 (1992). 

[16] Gh. Berger et. al, Phys. Lett. B227, 489 (1989). 

[17] M. Agiletta et. al, Europhys. Lett. 8, 611 (1989). 


20 












[18] M. Goodman et. aL, Nucl. Phys. (Proc. Suppl.) B38, 337 (1995). 

[19] The usage of the ‘ratio of ratios’ R as a valid indicator of the neutrino anomaly 
has been recently critically examined, see G.L. Fogli and E. Lisi, Institute for 
Advanced Study report lASSNS-AST 95/21. 

[20] H. Minakata, Preprint No. TMUP-HEL-9502, March 15, 1995; Preprint No. 
TMUP-HEL-9503, |hep-ph/9504222| , (1995). 

[21] L. Bento and J.W.F. Valle, Phys. Lett. B264, 373 (1991); J. Peltoniemi, A.Yu. 
Smirnov and J.W.F. Valle, Phys. Lett. B286, 321 (1992); J. Peltoniemi et. al, 
Phys. Lett. B298, 383 (1993); D.O. Galdwell and R.N. Mohapatra, Phys. Rev. D 
48, 3259 (1993). 


[22] See for example E. Ma and J. Pantaleone, Preprint No. UGRHEP-T140, March 
(1995); E. Ma and P. Roy, Preprint No. UGRHEP-T145, TIFR/TH/95-17 April 
(1995), E. J. Ghun, Anjan S. Joshipura and A. Yu. Smirnov, Report-no: IG/95/76, 


PRL-TH/95-7, |hep-ph/9505^ (1995). 

[23] For a review see L. Oberauer and F. Von. Feilitzsch, Rep. Prog. Phys. 55, 1093 
(1992). 

[24] B. Achkar et aL, Nucl. Phys. B434, 503 (1995). 

[25] F. Dydak et al, Phys. Lett. B314, 281 (1984). 

[26] B. Armbruster et al, Nucl. Phys. (Proc. Suppl.) B38, 235 (1995). 

[27] L. Borodovsky et al, Phys. Rev. Lett. 68, 274 (1992). 

[28] N. Ushida et al, Phys. Rev. Lett. 57, 2898 (1986). 

[29] M. de Jong et al., GERN-PPE/93-131 (1993); N. Armenise et al., GERN- 
SPSG/90-42 (1990). 


[30] P. Astier et al., GERN-SPSLG/91-21, GERN-SPSLG/91-48, GERN-SPSLG/P261 
Add.l (1991). 


[31] D. Saltzberg, hep-ph-9504343. 

[32] B.T. Gleveland et. al, Nucl. Phys. (Proc. Suppl.) B38, 47 (1995). 

[33] The theoretical rates quoted are from the last one in ref. [|], which is the latest 
in the series. 


[34] P. Anselman et. al, Phys. Lett. B285, 376 (1992); B327, 377 (1994). 


21 







[35] A.I. Abazov et. a/., Phys. Rev. Lett. 67, 3332 (1991); G. Nico et. a/., in Proc. of 
the XXVII Int. Conf. on High Energy Physics (Glasgow), Eds. P.J. Bussey and 
LG Knowles, Inst, of Phys. Publising, Bristol, p965 (1995); 

[36] K.S. Hirata et. al, Phys. Rev D44, 2241 (1991); Y. Suzuki et. al, Nucl. Phys. 
(Proc. Suppl.) B38, 54 (1995). 

[37] G.T. Ewan et. al, Sudbury Neutrino Observatory proposal. Pub. No. SNO-87-12 
(1987); Y. Totsuka, Tokyo Univ. Preprint No. IGRR-227-90-20 (1990), G. Arpe- 
sella et.al., Borexino Proposal, Eds. G. Bellini et. al., Univ. of. Milano (1991). 

[38] M. Kobayashi and K. Maskawa, Progr. Theor. Phys. 49 652 (1973); L. Maiani, 
Phys. Lett. B62, 183 (1976); U. Turke, et. al, Nucl. Phys. B258, 313 (1985); H. 
Harari and M. Leurer, Phys. Lett. B 181, 123 (1986). 

[39] E.D. Garlson, Phys. Rev. D34, 1454 (1986). 

[40] T.K. Kuo and J. Pantaleone, Rev. Mod. Phys. 61, 937 (1989) and references 
therein; See also D. Harley, T.K. Kuo and J. Pantaleone, Phys. Rev. D47, 4059 
(1993); M. Narayan et. al. Preprint No. IMSc/95-96/001 (1995). 

[41] A.S. Joshipura and P.I. Krastev, Phys. Rev. D50, 3484 (1994). 

[42] G.L. Fogli, E. Lisi and D. Montanio, Phys. Rev. D49, 3626 (1994). 

[43] V. Barger et al, Phys. Rev D43, 1110 (1991); Phys. Rev. Lett. 65, 3084 (1990); 
69, 3135 (1992); A. Acker, S. Pakvasa and J. Pantaleone, Phys. Rev. D43, 1754 
(1991); P.I. Krastev and S.T. Petcov, Phys. Lett. B285, 85 (1992); 299, 99 (1993). 

[44] V. Barger et. al, Phys. Rev. D43, R1759 (1991). 

[45] V. Barger, et. al, Phys. Lett. B93, 195 (1980); J. Phys. G 6, L165 (1980); V. 
Barger, K. Whisnant and R.J.N. Phillips, Phys. Rev D22, 1636 (1980); A. De 
Rujula et. al, Nucl. Phys.B168, 54 (1980). 

[46] J.J. Gomez-Gadenas and M.G. Gonzalez-Garcia, Preprint No. GERN-TH/95-80, 
April (1995). 


[47] S. Goswami, K. Kar and A. Raychaudhuri, Preprint No. GUPP/95-3, |hep- 
I ph/9505'^ (1995). 

[48] V. Barger and K. Whisnant, Phys. Lett. B209, 365 (1988). 

[49] A. Acker, J.G. Learned, S. Pakvasa and T.J. Weiler, Phys. Lett. B298, 149 (1993). 

[50] A. Acker, A.B. Balantekin and F. Loreti, Preprint No. Mad/NT/93-07, 
MAD/PH/# 774, July (1993). 


22 





[51] G. Barr, T.K. Gaisser and T. Stanev, Phys. Rev D39 3532 (1989); T.K. Gaisser, 
T. Stanev and G. Barr, D38, 85 (1988). 


23 




ii (a) ii (b) 


Fig. 1 


24 



